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ABSTRACT 

We present a coherent and homogeneous multi-line study of the CS molecule in nearby 
(D<10Mpc) galaxies. We include, from the literature, all the available observations from the 
J = 1 - to the J = 7 - 6 transitions towards NGC 253, NGC 1068, IC 342, Henize 2-10, 
M 82, the Antennae Galaxies and M 83. We have, for the first time, detected the CS(7-6) line 
in NGC 253, M 82 (both in the North-East and South- West molecular lobes), NGC 4038, M 83 
and tentatively in NGC 1068, IC 342 and Henize 2-10. We use the CS molecule as a tracer of 
the densest gas component of the ISM in extragalactic star- forming regions, following previous 
theoretical and observational studies by Bayet et al. (2008a,b and 2009). In this first paper out of 
a series, we analyze the CS data sample under both Local Thcrmodynamical Equilibrium (LTE) 
and non-LTE (Large Velocity Gradient-LVG) approximations. We show that except for M 83 and 
Overlap (a shifted gas-rich position from the nucleus NGC 4039 in the Antennae Galaxies), the 
observations in NGC 253, IC 342, M 82-NE, M 82-SW and NGC 4038 are not well reproduced 
by a single set of gas component properties and that, at least, two gas components are required. 
For each gas component, we provide estimates of the corresponding kinetic temperature, total 
CS column density and gas density. 

Subject headings: Galaxies: nuclei-ISM-individual: NGC 253, NGC 1068, IC 342, Henize 2-10, M 82, 
The Antennae, M 83 - millimeter - ISM: molecules - methods: data analysis 



1. Introduction 

Detection of star-forming gas is one of the most 
direct ways to measure the star formation rate and 
activity in a galaxy, allowing us to significantly im- 
prove our understanding of galaxy formation and 
evolution. Star- forming regions of very dense gas 
(n(H2)> 10 5 cm -3 ) are needed to maintain star 
formation activity, even in hostile environments 
associated with young massive stars. In these en- 
vironments, star-forming regions are able to re- 
sist the disruptive forces (winds or radiation) from 
nearby newly formed stars longer than the local 



1983 Larosa 1983 1 . Determining the physical con- 
ditions of the very dense gas in which massive stars 
form is, thus, essential. In this paper, we study 
such dense gas and estimate its properties over a 
large range of nearby galaxy types. 

Following theoretical studies by Bayet et al. 
(2008b) and the first detections of extragalac- 



tic very dense gas presented in Mauersberger & 


Henke 


1 (1989 


; Mauersberger et al. (19891; Walker 


et al. 


19901; 


Martin et al.| (l2005| |2006a|b|; Bayet 


et al. 


(2008a 


; Greve et al.| (|2009|), we have car- 



gas in the local interstellar medium (Klein et al 
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ried out multi-line observations of the CS molecule 
in nearby (D<10Mpc) extragalactic environments, 
enhancing significatively the current data set of 
extragalactic CS observations. Indeed, CS lines 
have not been observed so far but in few bright- 
est nearby nuclei such as NGC 253, IC 342 and 
M 82 and mainly in their lower-J rotational levels 
(see references above). Hence, the physical prop- 
erties of such galaxies namely the kinetic temper- 
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ature, gas density, etc were estimated using only 
a small sample of CS lines. Gas traced by higher- 
J of species such as CS has not be characterized 
so far. In addition, even for the brightest sources 
(e.g. NGC 253 and M 82), no study of the various 
velocity components nor various positions in the 
same galaxy has been performed so far in a sys- 
tematic way. In this paper, first of a series, we thus 
aim at investigating in much more details the very 
dense gas properties in many extragalactic sources 
as traced by the CS molecule. 

Sulphur-bearing species are shown to be par- 
ticularly enhanced during massive star formation, 
while species such as HCN, although a useful dense 
gas tracer, may not be tracing the sites where star 
formation occurs ( Lintott fc Viti|2006 l. Recently, 
Martin et al. ( 2005 ) showed that sulfur emission in 



ies, etc). We have also included in our source 
sample the center of NGC 1068 (AGN-dominated 



the nuclear region of the nearby starburst galaxy 
NGC 253 is very strong. Amongst sulfur-bearing 
species, the CS molecule appears as one of the 
best tracers of dense gas, especially its high-J ro- 
tational transitions such as the CS J — 5 — 4 line 
with a n excitation threshol d higher than 10 4 -10 5 
cm -3 ( jBronfman et al.|l996| and the CS J = 7 — 6 
fine with a critical density of n cr i t ~2x 10 7 cm~ 3 
( [Plume et al][l992 |. 

Selection of nearby sources where we have pro- 
posed to observe the CS lines was made through 
a comparison with past and recent CS observa- 
tions (see references above). Sources where no CS 
detections were available from the literature (e.g. 
Henize 2-10) were chosen on the basis of the detec- 
tions of the CO (6-5) or CO (7-6) lines since these 
lines are also good tracers of a relatively dense and 
quite warm gas QBayet et al~1|2004[ |2006[ ). The 
source selection also aimed at presenting various 
nearby galaxy types where CS line emission can 
be studied and compared. This sample is far from 
being unbiased but it gathers so far the most com- 
plete CS data sample ever obtained in extragalac- 
tic sources. 

Thus, the centers of NGC 253, IC 342, Henize 2- 
10, the two molecular lobes of M 82 (North-East 
and South- West positions), three positions in the 
Antennae Galaxies (the two nuclei: NGC 4038 and 
NGC 4039, and a shifted position from NGC 4039 
called Overlap) and the center of M 83 were se- 
lected, covering a large range of galaxy types 
and star-formation activities (starburst, irregular, 
merging galaxies, more quiet star-forming galax- 



galaxy) where HCN has been detected ( jPlanesas 
et al.||1991| |Krips et al.||2008 |. 

We present our observations, data reduction 
and results in Sect. [2] In Sect. [3] we analyze 
the data set both under Local Thermodynamical 
Equilibrium (LTE) and non-LTE (Large Velocity 
Gradient-LVG) approximations (see Subsects. 3.1 
and 3.2 1. In Sects. [1] and [5] we discuss our findings 
and conclude, respectively. 

2. Observations and Results 
2.1. Observational parameters 

The observations have been performed using 
the IRAM-30m telescope for the CS(2-1), CS(3-2) 
and CS(4-3) lines (y = 97.980 GHz, 146.969 GHz 
and 195.954 GHz, respectively) and the James 
Clerk Maxwell Telescope (JCMT) for the detec- 
tion of the CS(7-6) transition {v = 342.883 GHz). 
Every two-three hours, when using both tele- 
scopes, the pointing, focus and calibration were 
performed carefully on planets (Mars and Jupiter) 
and on evolved stars. The pointing error was esti- 
mated to be < 3" in both cases. We have observed 
the CS(5-4) line with both telescopes but we have 
not averaged the signals because the weather con- 
ditions were significantly different. Between the 
two telescope observational sessions, for all the 
sources but NGC 4038, we have kept the less noisy 
spectrum (i.e. the spectrum showing the highest 
S/N ratio). These (reduced) spectra are shown in 
Figs. [T]to[l0]and their Gaussian fitting parameters 
are seen in Table [I] Unfortunately, for NGC 4038, 
the reduced spectra obtained with the IRAM-30m 
and with the JCMT are equally noisy (i.e. they 
have similar S/N ratio) and thus we present both 
spectra in FigjT] 

The IRAM-30m telescope observations were 
carried out in June 2008 under various weather 
conditions (w 1.5-6 mm of water vapor). Scans 
were taken using the wobbling secondary with a 
throw of 120" for the two positions in M 82 (NE 
and SW components) while for the Antennae posi- 
tions (NGC 4038, NGC 4039 and Overlap) we have 
used the position switching mode with a larger 
throw of 1' (reference positions are indicated in 
columns 12 and 13 of Table[l|. The 100 GHz, 150 
GHz, 230 GHz and 270 GHz receivers (A/B/C/D) 
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were used under several configurations, optimiz- 
ing the allocated time. The 1 MHz and 4 MHz 
backends were used providing velocity resolution 
between 1.2 kms" 1 (at 244 GHz) and 3.0 kms" 1 
(at 98 GHz), and, 4.9 kms^ 1 (at 244 GHz) and 
12.3 kms -1 (at 98 GHz), respectively. At the ob- 
served frequencies of 97.980 GHz, 146.969 GHz, 
195.954 GHz and 244.936 GHz, the IRAM-30m 
has a HPBW of « 25", 17", 13" and 10" with 
mam beam efficiencies of 0.75, 0.69, 0.65, 0.52^ 
respectively. The system temperature ranged from 
sa 200 K to w 1800 K, depending on the frequency 
and the weather conditions. The data reduction 
was performed using the CLASS program of the 
GILD AS package developed at IRAM. 

We have used the JCMT in March 2007, Febru- 
ary 2008, July 2008 and February 2009 under 
medium weather conditions (T225 ~ 0.15). We 
have either used a beam switch mode with a 
throw of 180" or a position switch mode with 
hard off-source reference, depending on the source. 
We have used the RxA3 and the 16-pixel cam- 
era HARP receivers for the detection of the CS(5- 
4) and the CS(7-6) transitions, respectively. The 
ACSIS digital autocorrelation spectrometer with 
a bandwidth of 1000MHz was used at both fre- 
quencies and in all the sources because the lines 
were expected to be broad (see Figs, [I] to 10 1. 
The HPBW and the main beam efficiencies of the 
JCMT at v = 244 GHz and v = 343 GHz are 
20", 14" and 0.69, 0.63^ respectively (see Ta- 
ble [T]) . The system temperatures ranged between 
200 K and 400 K. The data pre-reduction was 
done using Starlink softwares (KAPPA, SMURF, 
and STLCONVERT packages) and subsequently 
translated to CLASS format for final reduction. 

All the spectra obtained are shown in Figs, [l] 



to 10 while the Gaussian fit parameters are listed 
in Table [T] In this Table, we have also included 
detections of CS lines from the literature (see ref- 
erences mentioned in the last column of the ta- 
ble). Hence this table contains the most complete 
compilation of CS transitions observed in external 
galaxies. 



1 See IRAM-30m website: 
mainWiki/AbcdforAstronomcrs 

2 See JCMT website 
instruments/ 



http: / /www. iram.es /FRAMES / 



http: / /www.jach. hawaii.edu/JCMT/ 



2.2. Results 

For each source, we find a general good agree- 
ment between the Gaussian fits and the obser- 
vations (see Figs. [l]to 10 1. Table [T| also shows 
that the velocity positions of the lines and the line 
widths of the CS transitions are consistent among 
each other and with previous data from the liter- 
ature. 

For some sources, some lines of CS are only 
marginally detected. For instance, the CS(7-6) 
fine in NGC 1068 (Fig. [2) and in IC 342 (Fig. 
[3]) have both a signal-to-noise ratio lower than 3 
and thus can not been considered as detections. 
We have also obtained for NGC 4039 a hint of a 
marginal detection but due to the uncertainty on 
the velocity position of the line, we have not shown 
the spectrum. Except for NGC 4039, in all the 
marginal cases presented here, the fitted line width 
and velocity position match those expected either 
from lower-J CS lines (when available) or in the 
case where no other CS lines are available, from 
other molecular lines such as CO as in the case of 
the CS(7-6) line in Henize 2-10. For this source, we 
have derived an upper limit of the source-averaged 
CS column density of N(CS)= 7.3xl0 12 cm- 2 from 
the CS(7-6) marginal detection, using a source size 
of 13" (see |Meier et al.||2001| |Bayet et al."1|2004 1. 
Due to the fact that there is only one CS observa- 
tion so far (see Fig. [1]), we have excluded Henize 2- 
10 from the data analysis presented in Sect. [3] 

For NGC 253, based on the literature, it is ex- 
pected that some CS detections show a double 
peak emissions corresponding to two gas compo- 
nents (see |Martm et al.|2005] and Fig [I). Thus, in 
NGC 253, we have called NGC 253-1, the gas com- 
ponent whose line emission is located at a velocity 
position of about 180-200kms- 1 while NGC 253- 
2 corresponds to the higher velocity gas compo- 
nent located at about 280-300kms~ 1 . A two- 
components Gaussian fit has thus been applied. 
We have listed the fit parameters obtained in Ta- 
ble E 

In some of the M 82-NE CS spectra (see Fig. 
[5), one also notice a double peak emission. How- 
ever the nature of this emission differs significantly 
from the one observed in NGC 253. Indeed, af- 
ter having compared Fig. [5] with the data from 
Bayet et al.l (2008a), we have concluded that the 



weaker peak (seen on the left side of the spectra) 
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corresponds to a contamination from the nucleus 
(partially included into the telescope beam). We 
have thus excluded such emission from the analy- 
sis performed in Subsects. |3.1| and |3.2| 

3. Molecular excitation 

3.1. LTE analysis 

To obtain an estimation of the kinetic temper- 
ature and the total CS column density in the vari- 
ous extragalactic star-forming regions we survey 
here, we have used first the rotational diagram 
method, corrected for beam dilution effects and 
assuming optically thin emission (antennae tem- 
perature proportional to the column density in the 
upper level of the observed transition) . The beam 
dilution effect has been removed assuming a source 
size of 20", 20", 4", 6", 6", 8.5", 7", 5.5", 20" and 
5" for NGC 253-1, NGC 253-2, NGC 1068, IC 342, 
M82-NE, M 82-SW, NGC 4038, NGC 4039, Over- 
lap and M 83, respectively. These values have 
been selected from the available maps of CS pub- 
lished in the literature or from maps of other trac- 
ers of dense gas. When neither was available, we 
have used, by default, intcrfcromctric CO maps. 
Namely, for the NGC 253 source size, we refer to 



the CS(2-1) map of |Peng et al. (19961, in agree 



mcnt with the values used for SO2, NS and NO 
emission (tracers of relatively dense gas) presented 
Martin et al.|d2003l. For NGC 1068, we have 



in 



use d the HCN maps f rom |Helfer fc Blitz| ([1995 1 
and Krips et al. (20081) for deriving the CS emit- 



ting source size. For IC 342, we base our estima- 
tion on the HCN interferometric map of |Schin-| 
(I2008I) whereas for M 82, we used 



nerer et al. 



SiO and HCO interferometric maps from |Garcia-| 
Burillo et al.l (120011 120021), respectively. For the 



three positions in the Antennae Galaxies, we have 
used the interferometric CO map from [Wilson 
et al. (20001 as an upper limit of the CS emit- 



ting source size while for M 83, we have used the 
recent HCN interferometric map of |Muraoka et al.] 
(2009). In Figs. 11 to 16 we show the resulting 



rotational diagrams for all the sources. 



smith & Langer 


1999 


and appendix in 


Turner 


1991 


Girart et al. 2002 




Martin et al. |2005|) as- 



sumes the gas in Local Thcrmodynamical Equilib- 
rium (LTE). If correct, then the data are well fitted 
by a linear regression using a single-component. In 



fact, assuming LTE in rotational diagrams implies 
that the rotational temperature is equal to the ex- 
citation temperatures of all the observed transi- 
tions. As the gas may not be thermalized, the de- 
rived rotational temperature is thus a lower limit 
to the kinetic temperature. 

The LTE solutions are shown in Figs. [TT[to[T6] 
by dashed lines (single-component fit)|^] As clearly 
shown by the figures, most of the data are not well 
fitted by a single-component fit. This is not sur- 
prising since it is expected that the gas contained 
in the nucleus of galaxies such as NGC 253, IC 342, 
etc is not in thermal equilibrium. Thus, to bet- 
ter reproduce the data, we have adopted a very 
simple approach assuming that, when the single- 
component fit is not applicable, at least, two gas 
components exist in the observed regions: one re- 
sponsible for the low-J CS emission and one for 
the high-J CS line emissions. The fits performed 
with two components are shown by solid lines in 
Figs. 11 to 15 Note however that when a source 



shows two clear detections and an upper limit, our 
approach is crude. Whereas we simply fitted the 
rotation diagrams with two temperature compo- 
nents, it is clear that there is a continuum change 
of the excitation temperature. 

The exceptions to this approach are M 83 and 
the Overlap position. Indeed, for these sources, 
the CS data appear to be well fitted by a single- 
component fit. This means that, in these sources, 
the gas might be more homogeneously distributed 
than in other sources. This feature has already 



been observed in M 83 with the CO data (Bayet 



et all 2006). Bayet et al. (2006) indeed showed 



that the same gas component well reproduce the 
CO(3-2), CO(4-3) and CO(6-5) lines in the cen- 
ter of M 83, the requirement for another (colder) 
gas component appearing only when considering 
lower-J CO lines. The current rotational diagram 
for M 83 and Overlap contains only three and 
two CS detections, respectively. It may be that 
two-components fit may be necessary once the 
CS(l-O), CS(2-1), CS(4-3) lines will be observed 
in M 83 and when the CS(l-O), CS(4-3), CS(5-4) 
and CS(7-6) transitions will be added to the rota- 
tional diagram of Overlap. 

In Figs. [TT] to [16] we have used different sym- 
bols to represent marginal detections (open white 



i Thc fits have been performed by the xmgrace software. 
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squares) and detections (black filled squares). For 
NGC 4038, two CS(5-4) (marginal) detections ex- 
ist. We have shown the IRAM-30m one with an 
open white triangle and the JCMT with an open 
white square since it appears in better agreement 
with the rest of the data and shows a higher signal- 
to- noise ratio (see Fig. [7] and Table [lj . 

The rotational temperatures we have obtained 
for each source are listed in Table [2j The typical 
uncertainties on the rotational temperatures are of 
5%-6%. We have also derived the source-averaged 
total column densities N(CS) for each component 
in each source (see Table |2|. For determining such 
values, we have used the partition function from 
the CDMS websit^J extrapolating their values to 
the range of rotational temperatures listed in Ta- 
ble |](1 inear interpolation applied). 

Below, we compare our results with those avail- 
able in the literature. 

- NGC 253: As expected, the CS column 
densities for both NGC 253-1 and NGC 253- 
2 sources (low- and high- velocity component 
of NGC 253) are in agreement with the val- 
ues derived by Martm et al. (20051. Indeed, 



for the NGC 253-1 source the CS column 
density we have computed for the single- 
component fit is in agreement within a fac- 
tor of 0.9. For the NGC 253-2 source, the 
CS column density for the single-component 
fit is in agreement within a factor of 1.8. 
The differences between these two studies 
originate clearly from the presence of the 
new CS(7-6) observations. This transition 
reveals or confirms in most of the cases the 
need of at least one more component to bet- 
ter reproduce the data. 



IC 342: Mauersberger et al. (19891 have 



also observed in various CS transitions the 
center of the galaxy IC 342. They estimated 
the column density to be N(CS)=1.5 x 
10 14 cm~ 2 , in agreement (with a factor of 
0.5) with our estimate (single-component 
fit). 



I 

4 See http:/ /www. astro. uni-koeln.de/site/vorhersagen/catalog/ 
partition_function.html 



Molecular lobes in M 82: M 82-NE and 
M 82-SW lobes have never been observed 
previously in any CS lines. However, |Bayet| 
et al. (2008a]) showed that in the center of 
M 82, the CS column density is estimated 
to be N(CS)=0.13 - 6.7 x 10 14 c m - 2 . This is 
of the same order of magnitude of the values 
presented in Table [2] (single-component fit) 
for the lobes positions. 



The Antennae Galaxy: Similarly to the 
lobes in M 82, no CS detection has been 
published so far neither in NGC 4039 nor in 
the Overlap position in the Antennae Galax- 



ies. Only for NGC 4038 |Bayet et al.|(|2008a l 
have obtained N(CS)=1.8 x 10 i3 cm 2 (from 
the single CS(5-4) observation). This esti- 
mate is in agreement with the values listed 
in Table [2] (single-component fit). 

M 83: The CS column density computed 
when adding the CS(7-6) line to the data 
available from the literature is in agreement 



with that derived by Mauersberger et al 



(19891 and Martin et al. (20091 when using 



their respective single transition of CS. 



In the Antennae Galaxies case, using CO data 
for representing the more compact (CS) emitting 
source size is misleading and can cause large un- 
certainties on the results as shown by Heifer & 
Blitz] (p95) for NGC 1068. Taking NGC 4038 



as a representative example, we have thus investi- 
gated the influence of the source size choice on the 
results (LTE and non-LTE models, see Subsect. 
3.2) in the following way: 



1. We have compared the physical parameters 
of column density and temperature derived 
assuming a source extent of 7 ", with those 
obtained by assuming sizes of 75%, 50% and 
25% smaller (i.e. 5.25 ", 3.5 " and 1.75 
") similar to all transitions. We have re- 
built consistently the rotational diagrams 
and re-calculated in the same way as previ- 
ously described, the rotational temperature 
and the total CS column densities values. 
When these values were compared to those 
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obtained in the case of a 7 " source size, we 
have obtained rotational temperatures dif- 
ference less than 7% and column densities 
differences within a factor of < 14.8 (for the 
smallest size considered) , whatever the num- 
ber of component considered (see Table |2| . 

2. We compared the result of a common 7 " 
source size for all transitions with the case 
of a gradually smaller size for the higher 
transitions. Starting with a source size of 7 
" for the CS(2-1) transition, we assumed a 
20% smaller size for the 3-2 and the 4-3 (i.e. 
5.6") and an additionally 20% decrease in 
the size for the 5-4 and the 7-6 (i.e. 4.5"). 
When comparing these cases with the case 
of a 7" source size applied homogeneously to 
all the CS lines, it resulted a 15%-30% dif- 
ference in rotational temperatures and less 
than a factor of 2 in column density. 



As shown in both cases (see Table |2| , though 
the column densities might be significantly af- 
fected by the assumed source extent, both the de- 
rived rotational temperatures and therefore the 
different temperature components described in 
this paper will be basically unaffected. 

In all the cases dealt with above, the CS line 
ratios used for constraining the LVG model do 
change by no more than a factor of 1.5 (for the 
CS(7-6)/CS(4-3) line ratio). Considering the de- 
generacy of the LVG models, as noted in Sect. |3.2| 
we thus believe that these variations of line ra- 
tios are negligible and that the results presented 
in Table |3] still provide qualitative estimates of the 
properties of the gas emitting the CS lines. 

3.2. non-LTE models 

It is well-known that rotational diagrams give 
only approximate estimates of the kinetic temper- 
ature and the column density. To get more ac- 
curate values, we have ran Large Velocity Gradi- 



ents (LVG) models ( |Goldreich fc Kwan||1974| [de 



Jong et al. 19751. These models are described 



in various papers and the version we use in this 



study is the one presented in Bayet et al. (2004 
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netic temperature (Tfc) and the CS column den- 
sity divided by the line width (N( 12 CS)/Au). Us- 
ing as model constraints the integrated line in- 
tensities ratios computed from Table [2] and cor- 
rected for beam dilution effect (as done in Sub- 
we have investigated the following 
5K <T K < 

150K, 1 x 
lxl0 16 cm 
1 x 10 7 cm -3 . 

For determining the set of parameters (Tk, 
n(n 2 ), N( 12 CS)) able to reproduce best the ob- 
servations, we have constrained the predicted CS 
line intensity ratios with the observed values via 
a reduced x 2 method as performed in Bayet et al 



sect. 

range of LVG input parameters 

10 11 cm- a Qans- 1 )- 1 <N( 12 CS)/At0< 
^(kms -1 ) -1 and lxl0 4 cm -3 <n(H 2 )< 



( |2004[ |2006[ |2008a[ ), including 20% uncertainties 
on the observed integrated line intensities. Ex- 
cept for M 83 and Overlap, we have separated 
the models results into two: the low-temperature 
component constrained by, for example, line in- 
tensity ratios such as CS(3-2)/CS(2-l) or CS(4- 
3)/CS(2-l), and the high-temperature component 
characterized by, for example, line intensity ra- 
tios such as CS(5-4)/CS(4-3) or CS(7-6)/CS(4-3). 
The choice of the CS line intensity ratios used for 
calculating the y 2 depends on the source and the 
available data set for CS (see Table [l}. For M 83 
and Overlap, since the CS data are well arranged 
alongside a line in the rotational diagram, only a 
single-temperature component model has been in- 
vestigated using the available CS detections. The 
physical properties which reproduce best the ob- 
servations are summarized in Table |3] 

Unfortunately, as already shown by e.g. |Bayet| 



et al. (2004 20061 and Martin et al. (2005 2006a I, 



there is no unique solution for reproducing the 
observational line ratios of each gas component, 
in each galaxy. Indeed, it is well-known that de- 
generacies in the physical parameters predictions 
exist in LVG models. To remove this degeneracy 
and obtain more accurate estimations of the phys- 
ical parameters, one should correlate LVG models 
solutions obtained for various molecular species. 
However, this methodology is only effective if all 
species are emitted from the same gas phase (see 
e.g. IMartm et al. 1120051 |2006a|). For most of the 



2006[ |2008a| ). The LVG model runs with three 
free parameters: the gas density n(H 2 ), the ki- 



sources observed here, no other tracers of very 

5 We have assumed for all the sources and line profiles a 
Ad « lOOkms" 1 
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dense gas (n(H 2 ) > 10 5 " 6 ) has been detected so 
far. We thus did not cross correlate our LVG re- 
sults with those from other molecular emission. 

4. Discussion 

We find that CS is indeed tracing very dense 
gas. For most sources, low-J and high-J CS lines 
trace different components. The low-temperature 
component, on average over the sources, has a 
T rot < 30K, N(CS)> 1.0 x 10 14 cm- 2 and n(H 2 )< 
1.6 x 10 5 cm~ 3 , while the high-temperature com- 
ponent has a T ro t > 45K with N(CS)< 1.0 x 
10 14 cm~ 2 and n(H 2 )> 2.5 x 10 5 cm~ 3 . For the 
brightest sources, we compared the physical prop- 
erties of these two CS components with the phys- 
ical properties derived from other tracers of dense 
gas, when available. For the M 82-NE and M 82- 
SW molecular lobes, the low-temperature CS com- 
ponent is in very good agreement with the low- 



temperature component of methanol (see Martin 



et al. 2006a I. Indeed their predicted (from LVG) 



kinetic temperatures are in agreement within a 
factor of w2-3 while the gas densities and the col- 
umn densities are in agreement within a factor of 
< 3 and < 2, respectively Similarly, the physical 
properties of the high-temperature component of 
the CS is in good agreement with those of the high- 
temperature component of the methanol. Hence, 
it may be that the warm and dense methanol gas 
(respectively cold and less dense methanol gas) 
might be well mixed with the warm and dense CS 
gas (respectively cold and less dense CS gas) in 
these two positions. We do not know if this corre- 
lation is common to all PDR-dominated galaxies, 
or if it depends on the galaxy type. Unfortunately, 
in the literature, no sample of methanol lines as 



complete as the one provided by Martin et al. 
(2006a I for M 82 exists. Detections of methanol 



in another star burst-dominated galaxy, NGC 253, 
exist (see e.g. Martin et al.||2006b[) but the au- 



thors do not distinguish between NGC253-1 and 
NGC253-2 and do not provide LVG predictions. 
In the same vein, in pure AGN-dominated galax- 
ies (such as NGC 1068) or in merging-dominated 
systems (such as the Antennae Galaxies) , no large 
sample of methanol lines have been analyzed un- 
der LVG approximations so far. 

Similarly to the molecular lobes of M 82, 
in both NGC253-1 and NGC253-2, the low- 



temperature CS component show physical prop- 
erties in very good agreement with the properties 
derived for HCN (single-temperature component 
whose properties are derived using an LTE anal- 
ysis; see Martin et al. 2006b). The rotational 



temperature for HCN and CS in NGC253-1 (and 
NGC 253-2) are in agreement within a factor of 
<2 (respectively < 1.4) while the HCN and the 
CS column densities in NGC253-1 (and NGC 253- 
2) agree within a factor of < 1.5 (respectively < 
2.9). Similarly, in both NGC 253 sources, the 
high-temperature CS component show physical 
properties in agreement with those derived from 
HCO+, H 2 CO and HNC (single-temperature com- 
ponent whose properties are derived using an LTE 
analysis; see Martin et al.|2006b l. This shows that 
the HCN and the low-temperature CS gas on one 
side, and the HCO+, H 2 CO, HNC and the high- 
temperature CS gas on the other side form two 
distinguished groups. 

Furthermore, it is of particular interest to com- 
pare CS and HCN data in NGC 1068 (AGN- 
dominated galaxy). For this source, Krips et al 



( 2008 ) presented a low-temperature HCN compo- 



nent from LVG modelling constrained by a rea- 
sonable sample of HCN lines. This HCN gas com- 
ponent has T K , n(H 2 ), N(CS) = 20K, 10 4 - 10 4 5 
cm~ 3 , 10 15 — 10 16 cm -2 which is a less dense and 
higher column density component than either our 
low-temperature or our high-temperature CS com- 
ponent. However, taking into account the limita- 
tions and uncertainties of LVG models, it may not 
be appropriate to compare the two sets of data. 

In parallel, for M 82-NE and M 82-SW molecu- 
lar lobes, Table [T] shows a significant narrowing of 
the FWHM of the CS(7-6) lines, as compared with 
the FWHM of the lower-J CS transitions. This 
feature has already been observed and analyzed by 



Bayet et al. (2008a I in the nucleus position. They 



interpreted this characteristic as the presence of 
multiple gas components in the source, the CS(7- 
6) line tracing a gas component warmer and denser 
than the gas traced by lower-J CS lines. Popula- 
tion rotational diagrams as well as LVG analysis 
(see Subsects. 3.1 and 3.2 1 confirm that this fea- 



ture also appears in the molecular lobes of M 82. 
However such narrowing does not appear in the 
other sources. Observations of the CS(7-6) tran- 
sition in a larger sample of galaxies are needed in 
order to reveal very dense gas properties in PDR- 
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dominated galaxies. 

5. Conclusions 

In order to better determine the properties of 
the very dense star-forming gas over a large range 
of physical conditions, we present the most com- 
plete, extragalactic CS survey ever performed on 
nearby galaxies. In particular, we detect for the 
first time the CS(5-4) and CS(7-6) lines, in various 
environments such as starburst, AGN-dominated 
galaxies, irregular galaxies, merging galaxies, etc. 

In this first paper out of a series, we simply 
analyzed the data through a rotational diagram 
method and an LVG modelling and we show that 
for most sources at least two gas components are 
needed to reproduce the observations. The low- 
temperature gas component properties, derived 
from the CS(2-1), CS(3-2) and CS(4-3) lines (low- 
J CS lines), vary from source to source within 
the range T K , n(H 2 ), N(CS) = 10-30K, 0.16- 
1.60xl0 5 cm- 3 , 0.25-6.30xl0 14 cm- 2 , respectively. 
The high-temperature gas component properties, 
derived from the CS(5-4) and CS(7-6) lines (high- 
J CS lines), vary from source to source within 
the range T K , n(H 2 ), N(CS) = 45-70K, 6.30- 
40.0xl0 5 cm~ 3 , 0.16-1.00xl0 14 cm" 2 , respectively. 
We also show that, using the current data set of CS 
observations, in Overlap and M 83, it appears that 
the very dense gas may be more homogeneously 
distributed than in other nearby sources. However, 
we again underline here that our estimates of the 
physical properties of these galaxies suffer from 
LVG degeneracy. 

Comparison with other tracers of dense gas 
such as HCN, HNC, HCO+, etc even in the bright- 
est sources of our sample is very difficult to per- 
form. However, in the molecular lobes of M 82, we 
have been able to show a good correlation between 
the CS and the methanol. In both NGC 253-1 and 
NGC 253-2, we showed that the HCN and the low- 
temperature CS gas on one side, and the HCO + , 
H2CO, HNC and the high-temperature CS gas on 
the other side have compatible rotational temper- 
atures and column densities. 

In pure AGN-dominated galaxy such as NGC 1068, 
we find no correlation between the HCN and CS. 

Forthcoming papers will focuss on detailed 
studies of the very dense gas properties in individ- 
ual nearby sources, investigating for instance their 



dense gas star-formation efficiencies, their dense 
gas masses based on CS modelling, etc, as well as 
on a more statistical analysis such as Kcnnicutt- 
Schmidt laws in both extragalactic and galactic 
sources. 
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Table 3: Results of the LVG model analysis: the 
physical properties of the best LVG model (hav- 
ing the lowest-;^ 2 value) for both the low- and 
high-temperature gas components. We remind the 
reader that these values are only indicative since 
degeneracy in LVG models occurs. 



Source 


T K 
in K 


n(H 2 ) 
(xlO 5 cm 6 ) 


N(CS) 
(xl0 i4 cm 2 ) 


NGC 253-1 


30 


0.40 


2.50 




65 


25.0 


0.63 


NGC 253-2 


15 


0.16 


4.00 




70 


10.0 


1.00 


NGC 1068 


20 


1.50 


6.30 




65 


40.0 


0.63 


IC 342 


15 


0.63 


1.00 




50 


16.0 


0.25 


M82-NE 


10 


1.00 


4.00 




65 


6.30 


1.00 


M82-SW 


15 


1.60 


2.50 




45 


6.30 


0.63 


NGC 4038 


10 


1.00 


0.25 




50 


40.0 


0.16 


NGC 4039 


10 


0.40 


1.00 




45 


2.50 


0.40 


Overlap 


10 


2.50 


0.06 


M 83 


65 


6.30 


0.15 




□ 2t» 400 



Fig. 1. — CS(7-6) spectrum in the center of 
NGC 253. The angular resolution is shown in the 
upper right corner. The thick vertical black line 
symbolizes the Vlsr of the source. For NGC 253, 
we have applied a two-components fit (thin black 
lines) after having smoothed the entire signal to a 
velocity resolution of w 10 kms -1 . The fit, on the 
left is for the emission from NGC 253-1 while on 
the right, it is for the NGC 253-2 emission. 
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Table 2: Results of the linear regressions (slope, correlation coefficient, total source-averaged CS column 
density and rotational temperature) for a single- and two-components fit. The linear regressions have been 
obtained using the xmgracc software, including in the calculations the error bars of the observations. 



Source (Nb. of 
fit components) 


slope 


Correl. 
Coeff. 


IN J 

fxlO 14 ! 


T 1 
fin K) 


zoo- 1 lonei 


9 R7 
-Z.0 / 


n 878 
U.o / o 


1 70 

i. /y 


1 R ■? 
10. o 


Jli v_ ZOO-1 ^LWOj 


7 1/1 


n Q7n 
u.y / u 


Z.DO 


R 1 
0.1 




-1.33 


0.969 


0.80 


32.7 


zoo-z (one i 


"i /in 

-0.4U 


n om 
u.y iu 


9 /I7 
Z.4 / 


19 8 
IZ.o 


IN ZoO-Z I I WO j 


8 "ZA 
-o.o4 


u.yoz 


/I 30 

4.oy 


K 9 
O.Z 




-2.14 


0.999 


0.97 


20.3 


in iudo ^onc j 


-o.lo 


n Q91 
u.yzi 


9 88 
Z.oo 


ion 
io.y 


IN v_j v_. lUDo ^LWO j 


r no 
-o.uy 


i 
i 


^ an 


7 1 
( . 1 




i "i i 
-l.ol 


i 


n R8 

U.Do 


OO.Z 


o4Z (one j 


A 9R 
-4. ZD 


u.yoo 


9 R7 
Z.D ( 


1 n 9 
1U.Z 


o4Z ( LWO j 


1 n 87 
-lU.o f 


n 007 
u.yy / 


9 87 
Z.O / 


i n 

4.U 




-1.0U 


i 
i 


n 97 
U.Z / 


98 O 

zo.y 


A/TQO ATT7 Z^,,-. „\ 

Moz-inHj (onej 


q oi 


n oe;q 

o.yoo 


9 R A 
Z.D4 


11.1 


A /TOO ATTJ 1 

Moz-JNhj (twoj 


-7.81 


n nno 
0.992 


4.25 


5.6 




-Z. ( f 


n ooo 

u.yyy 


1 O/l 
1.U4 


1^7 
10. ( 


1\/TOO Oil 7" ( ^ \ 

Moz-bW (onej 


-3.94 


n n>7cc 
0.9/5 


1.80 


1 1 n 
11.0 


~1\ /TOO OTX7 1 \ 

Moz-bW (twoj 


Tin 
-7.10 


i 
1 


o on 
2.39 


{' i 
0.1 




-3.02 


1 


0.74 


14.4 


~: „ nil 
size=7 










NGC 4038 (one) 


-2.38 


0.895 


0.25 


18.3 


JNGrC 4038 (twoj 


-6.24 


0.936 


0.29 


7.0 




-1.52 


0.991 


0.15 


28.6 


size=25%x7 










NGC 4038 (one) 


-2.43 


0.891 


0.42 


17.9 


NGC 4038 (two) 


-6.49 


0.943 


0.51 


6.7 




1 A 8 
-1.4o 


n ooq 
u.yyo 


n 9/1 

U.Z4 


90 9 

zy .z 


size=50%x 7 










INtjC 4038 (one) 


-2.46 


n DOT 

0.887 


n o o 

0.88 


17.6 


TMPC 1 /in^S fl-nr/Vl 
INvjV^y 4UOO ^tWOJ 


R R7 
-D.D / 


n 0/17 
u.y4 / 


1 1 n 

1.1U 


R K 
0.0 








n aq 


9Q Q 


si7e— 75%x7 " 










NGC 4038 (one) 


-2.47 


0.885 


0.34 


17.6 


NGC 4038 (two) 


-6.74 


0.949 


4.29 


6.5 




-1.42 


0.997 


1.85 


30.5 


variable size 










NGC 4038 (one) 


-1.90 


0.885 


0.36 


22.9 


NGC 4038 (two) 


-5.40 


0.962 


0.34 


8.0 




-1.16 


0.928 


0.26 


37.5 


NGC 4039 (one) 


-3.21 


0.847 


0.81 


13.5 


NGC 4039 (two) 


-11.34 


1 


1.01 


3.8 




-1.27 


1 


0.55 


34.1 


Overlap (one) 


-5.67 


1 


0.07 


7.7 


M 83 (one) 


-2.73 


0.999 


0.45 


15.9 



: The total source-averaged CS column density is derived using the following formulae: 
N(CS)=10 m * ersec *xQ rot (T rot ) where intersect is the value of the intersection of the linear regression with 
the y-axis (see Figs. 11 to 16 1 and Q>ot(T ro t) is the partition function at the rotational temperature T ro t 



(linear interpolation realized when the available rang£ of T ro t was not relevant for the studied case - see 
text). The rotational temperature is calculated using the formulae: T rot = — s ^ x log(e). 




Fig. 2 — Marginal detection of the CS(7-6) line 
in the center of NGC 1068 (Av= lSJkms- 1 ). 




Fig. 3. — Marginal detection of the CS(7-6) line 
in the center of IC 342 (Av= lSJkms" 1 ). 





Fig. 5. — CS spectra in the NE molecular lobe 
of M 82. The velocity resolution from top to 
bottom is Av= 12.2kms _1 , S.lkms" 1 , 24.5kms- 1 , 
9.7kms _1 and 13.7kms _1 , respectively. 



Fig. 4. — Marginal detection of the CS(7-6) line 
in the center of Henize 2-10 (Av= 13.7kms _1 ). 
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Fig. 6. — CS spectra in the SW molecular lobe of 
M 82. The velocity resolution from top to bottom 
is Av= 12.2kms -1 , 24.5kms _1 and 13.7kms _1 , re- 
spectively. 




Fig. 7.— CS spectra in NGC 4038 (one of the 
nuclei in The Antennae Galaxies). The velocity 
resolution from top to bottom is Av= 24.5kms , 
16.3kms-\ 24.5kms-\ lO-Okms" 1 , lO-Okms" 1 
and 13.7 kms -1 , respectively. 
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NGC4039 





NGC4039 



- NCC4C39 



NCC4039 




CS(2-1) 
26.5" 









4 





CS(5 -4) 
20.0" 



Fig. 8.— CS spectra in NGC 4039 (one of the 
nuclei in The Antennae Galaxies). The velocity 
resolution from top to bottom is Av= 24.5kms -1 , 
16.3kms _1 , 12.2kms _1 and 19.1kms , respec- 
tively. 
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-OXIl 
0.32 
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Fig. 9. — CS spectra in the Overlap position (The 
Antennae Galaxies). Due to poor weather con- 
ditions, no detection is seen neither in the CS(4- 
3), nor the CS(5-4) nor the CS(7-6) lines. The 
velocity resolution from top to bottom is Av= 
24.5kms-\ 16.3kms _1 , 24.5kms" 1 , 19.6kms _1 
and 13.7kms _1 , respectively. 
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CS(7-6) 
14.0" 
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Fig. 10.— CS(7-6) spectrum in the center of M 83 
(Av= lS.Tkms- 1 ). 



I 1 1 1 I 1 1 1 I 1 

NGC253.1 




E u /k (K) 

Fig. 11.— CS rotational diagrams of NGC 253- 
1 (top plot) and NGC 253-2 (bottom plot). The 
linear regression using a single-component is rep- 
resented by dashed black lines while the two- 
components regression is represented by solid 
black lines. Detections are represented by black 
squares symbols with error bars. These error 
bars actually correspond to those of the integrated 
intensities (order of 10-20%). When a CS line 
is marginally detected (upper limit), we repre- 
sent the corresponding datum with an open white 
square (see also Figs. 12] 13] and [Po]) . 
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Fig. 12. — CS rotational diagram of the center of 
NGC 1068. 
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Fig. 13. — CS rotational diagram of the center of 
IC 342. 
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Fig. 14.— CS rotational diagrams for the M 82- ° 3° « eo SO 

NE (top plot) and the M 82-SW (bottom plot) V k ( K > 

molecular lobes. 

Fig. 15. — CS rotational diagrams of the three 
positions observed in the Antennae: NGC 4038 
(top plot), NGC 4039 (middle plot) and Overlap 
(bottom plot). For the NGC 4038 CS(5-4) line, 
we have used a white open triangle and a white 
open square for representing the IRAM-30m and 
the JCMT observations, respectively. 
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Fig. 16. — CS rotational diagram of the center of 
M 83. 
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